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Abstract 

The thermodynamics and critical exponents and ampUtudes of 
high temperature and dense matter near the chiral critical point are 
studied. The parameterized equation of state matches on to that cal- 
culated with lattice QCD at zero chemical potential and to the known 
properties of nuclear matter at zero temperature. The extent to which 
finite size effects wash out the phase separation near the critical point 
is determined. 
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1 Introduction 



The up and down quark masses are very small but not zero. In consequence, 
the conventional wisdom is that there is no true thermodynamic chiral phase 
transition at finite temperature T and zero baryon chemical potential fi. 
Instead, there is expected to be a curve of first-order phase transition in the 
/x-T plane that terminates in a second-order phase transition at some critical 
point {fic,Tc). The location of the critical point is obviously of quite some 
interest. This topic has been under intense theoretical study using various 
effective field theory models, such as the Namu Jona-Lasinio model [I]- [3], a 
composite operator model a random matrix model [5], a linear a model 
[3j, an effective potential model [H], and a hadronic bootstrap model [7], 
as well as various implementations of lattice QCD [8j-^llj. Reviews of this 
subject were presented by Stephanov [12] and Mohanty [13]. 

This subject is also of great interest because collisions between heavy 
nuclei at medium to high energy, such as at the future Facility for Antipro- 
ton and Ion Research (FAIR), or possible low energy runs at the Relativistic 
Heavy Ion Collider (RHIC), may provide experimental information on the 
phase diagram in the vicinity of a critical point. One characteristic signa- 
ture would be large fiuctuations in phase space of conserved quantities, such 
as charge, baryon or strangeness, on an event by event basis [T^ [T5] . The 
variance of the distributions is proportional to the spatial size of the corre- 
lated region, which could be rather small due to the finite size and lifetime 
in heavy ion collisions [16]. This led to the suggestion to measure higher 
moments to search for non-Gaussian behavior [17]. In order to study these 
effects quantitatively, not only are measurements needed but also dynamical 
simulations of phase separation and fluctuations in heavy ion collisions |18j . 

The goal of this paper is to understand the basic features of the equa- 
tion of state near the QCD chiral critical point and the magnitude of phase 
fluctuations in its vicinity. The essential requirements are to incorporate the 
critical exponents and amplitudes, and to match on to lattice QCD results 
at /X = and to nuclear matter at T = 0. We will accomplish this by param- 
eterizing the Helmholtz free energy as a function of temperature and baryon 
density so as to incorporate the above requirements. 

This work is similar to the study of the nuclear liquid gas phase transition 
in [19]. An obvious difference is that the latter studied a transition between 
nuclear liquid and gas, whereas the present paper will study the transition 
between quark matter and hadronic matter. Apart from that, this paper will 
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develop a description of the equation of state that has the correct critical 
exponents which are not integers or simple fractions; in other words, it is not 
a mean field theory. It will also include the critical amplitudes, which are 
universal. The chiral phase transition is in the same universality class as the 
3D Ising model and liquid gas phase transition. Finally, the parameterization 
will incorporate knowledge about the zero baryon density equation of state, 
computed by lattice QCD, as well as knowledge about the high density be- 
havior of nuclear matter at zero temperature. Perhaps the closest work that 
addressed these issues was ref. [20] which blended a parameterization of the 
3D Ising model equation of state into quark and hadron equations of state. 
In that work there was an ambiguity as to how to relate the magnetization 
to the baryon chemical potential. In this paper there is no such issue. These 
two parameterizations can perhaps be viewed as alternatives which provide 
some idea as to the range of uncertainty in how to describe matter near the 
chiral critical point. 

The outline of the paper is as follows. In Sec. 2 we summarize some 
thermodynamic relations and definitions of the relevant critical exponents. 
In Sec. 3 we develop a relatively general parameterization of the equation 
of state in the vicinity of the critical point whose motivation comes from 
mean field theory. In Sec. 4 we fix the parameters to best match onto what 
is known about the equation of state at finite temperature but zero baryon 
density, and at finite baryon density but zero temperature. In Sec. 5 we 
show the numerical results obtained with the developed parameterization 
and parameters. In Sec. 6 we apply the Landau theory of fluctuations away 
from the stable thermodynamic phases to estimate the magnitude of density 
fluctuations, which turn out to be surprisingly large. Concluding comments 
are made in Sec. 7. Those wishing to know only the results should read Sec. 
5 and the Appendix which summarizes the parameterized equation of state. 

2 Thermodynamic Relations 

Consider the equation of state of matter in the vicinity of a critical point 
which is in the same universality class as the liquid-gas phase transition 
and the 3D Ising model. It is advantageous to work with the Helmholtz free 
energy density / which is a function of the temperature T and baryon density 
n. Some useful thermodynamic relations involving the pressure P, baryon 
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chemical potential entropy density s and energy density e follow. 

/ = f{n,T) (1) 

. ^ (4) 

e = f{n,T) + Ts{n,T) (5) 

Clearly these satisfy the thermodynamic identity e = —P + Ts + nn. The 
heat capacity per unit volume cy and isothermal compressibility are 

cv - t'-^ (6) 

The baryon number susceptibility is xb — ri^i^T- In what follows we shall 
focus on kt rather than xb since they are so directly related to each other. 

When discussing a critical point with critical temperature Tc and baryon 
density ric it is useful to define 

T -T 

t - (8) 



The critical exponents a,/3,7,5 are defined as follows. When t ^ |?7| and 
t > 

Cv ~ t"" (10) 
Kt - r^. (11) 

Along the coexistence curve 

ni-Ugr. {-tf . (12) 

Along the critical isotherm 

P-Pc^ |7;|''sign(7;) . (13) 
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Mean field theories normally give a = 0,/3 = l/2,7 = l, and S = 3, as we 
shall see. Typical fluids are measured to have a ^ 1, /3 ^ 1/3, 1.2 < 7 < 1.3, 
and 4 < (5 < 5. For experimental measurements and results see [_21j and 
references therein. The 3D Ising model has a = 0.11, /3 = 0.325, 7 = 1.24, 
and 6 = 4.815. A good general reference is [22]. We will now proceed to 
parameterize the equation of state near the chiral critical point of QCD at 
increasing levels of sophistication. 

3 Parameterizing the Equation of State 

As mentioned earlier, there are certain properties near a critical point that 
are identical for all theories within the same universality class. Away from 
the critical point the equation of state depends on the details of the degrees 
of freedom and the interactions among them. To proceed we will flrst briefly 
review generic descriptions that arise from most if not all mean flelds theories. 
The results will be used to motivate a more sophisticated parameterization 
relevant for QCD. 

3.1 Mean field theories 

By mean fleld theories it is meant that although interactions are included, 
correlations among the particles are not. This usually results in thermody- 
namic variables scaling as rational powers of rj and t. This can be represented 
by expanding / in a Taylor series in rj about rj = 0. 

f = j:mv' (14) 

A:=0 

The coefficient functions fk{t) themselves may be expanded in a Taylor series 
about t = and they all have energy dimension 4. The resulting pressure is 

k=0 

Pk = {k + l)U+i{t) + {k-l)fk{t). (15) 

At the critical point both dP{n,T)/dn = and d'^P{n,T) /dn^ = 0. This 
implies that Pi(0) = and P2{0) = 0, or equivalently, /2(0) = and /3(0) = 
0. Similarly the other thermodynamic variables may be obtained, such as 

= -E(^ + 1)^+1 W^' (16) 

""-c k=0 
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and 

s--^i:mv' (17) 

where the prime denotes differentiation with respect to t. 

The simpfest model is to set fk{t) = for all A; > 4. A quartic polynomial 
for the free energy is typical in particle and condensed matter physics. This 
means that P is quartic and /i is cubic in r], which allows for the usual S- 
shaped curves. Phase coexistence requires that the pressures and chemical 
potentials of the two phases be equal for t < 0. 

Pir)i,t) = Pir)g,t) (18) 
KVl,t) = KVg,t) (19) 

The subscripts / and g stand for the liquid high density phase and the gaseous 
low density phase, respectively. This determines r]i{t) > and r]g{t) < along 
the coexistence curve. For this model 

2/2 ivi - %) + (3/3 + /2) {vf - Vl) + (4/4 + 2/3) {vf - rfg) + 3/4 {vt - v'g) = 

(20) 

and 

2/2 {r]i - %) + 3/3 {vf - V^) + 4/4 {vf - Vl) = (21) 

Since both /2 and vanish at the critical point, should not in order that 
Tji and rjg go to zero as i — > 0. These equations apparently do not have any 
simple solution. Therefore assume that fsit) = 0, which is certainly not the 
most general case but it does allow us to gain valuable insight. Then the 
solution to these equations is 



Vi{t) = -%{t) = 

So the function /2(t) should be negative for t < and positive for t > 0. 
The simplest choice, usually obtained in mean field approximations, is that 
f2{t) ~ t. This gives f3 = 1/2. The resulting coexistence curve in the T- 
n plane is symmetric about ric- Real fluids oftentimes have an asymmetric 
curve. 

Along the coexistence curve the chemical potential is 

ncfi.{T) ^ hit) (23) 



(22) 
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If, for example, the coexistence curve T versus // is parameterized as an 
ellipse then the function fi{t) is determined. The coexistence pressure is 

/'.(T) = P„W + || = -/„M+/:W + f| (24) 
Along the critical isotherm 

P{n, Tc) - PK, T,) = 4U{0)v' + 3h{0W (25) 
so that the critical exponent S — 3. The thermal compressibility is 

Pi(t) + 2P2it)ri + ?>P^{t)rf + 4P4(i)r7'] (26) 

When i > Ir^l ^ and i > 

so that if f2{t) vanishes linearly then the critical exponent 7=1. In the 
same limit 

cv = -^m) (28) 

c 

If foit) is a regular function then the critical exponent a — Q. 

The limit of meta-stability is the isothermal spinodal. It is determined 

by the condition dP{n,T)/dn = 0. In this model one finds, with f^it) = 0, 
that the lower limit is % = r]g/ < and the upper limit is 772 = Vi/VS > 0. 
For rjg < Tj < rji the system is in a meta-stable gas phase, and for rj2 < r] < rji 
the system is in a meta-stable liquid phase. 



3.2 A realistic parameterization 

Now we construct a model that has the correct critical exponents. The 
most important consideration is to obtain the correct value of 5 which is an 
irrational number. Motivated by the mean field theories, we parameterize 
the free energy as 

/ = m + h{t)r) + h{t)rf + UitM'' . (29) 
The pressure, chemical potential, and entropy density arc 

P = -/o + /i + 2/277 + hn^ + afM"'^ s^s^iv) + (^ - i)M^r (30) 
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n,i^ = /i + 2f2V + afM"'' sigii(77) (31) 

TcS^-r^-An-f^rf-fM"- (32) 

Phase coexistence is determined by equal pressures and chemical potentials 
at the same temperature but different densities. 

2h{vi-V9) + ^U[Wr' + W) = (33) 
f2{ri-rf^-r{a-l)UM^-\rj,Y) = (34) 

The second of these has an obvious solution for rig — —rii. When substituted 
in the first equation we find 



1 

0--2 



(35) 



Since /2(0) = the pressure along the critical isotherm is 

P(n,Tj - P(n„T,) = a UOM"'' sign{r]) + (a - l)/a(0)|r/r (36) 

and so the critical exponent 6 = a — 1. 

The limit of met a- stability is the isothermal spinodal. It is determined by 
the condition dP{n, T)/dn — 0, as mentioned earlier. Now one finds that the 
lower limit is rji — rjg/ 5^l^^~^'> < and the upper limit is 7)2 = r]i/5^'^^-^'^ > 0. 
For Tjg < Tj < r]i the system is in a meta-stablc gas phase, and for r]2 < rj < r)i 
the system is in a meta-stable liquid phase. In the range rji < rj < rj2 the 
system is unstable against isothermal fluctuations. 

The density difference goes to zero as 

{-tr ■ (37) 

In the 3D Ising model and in real liquid-gas transitions it turns out that the 
thermal compressibility kt diverges as K+t~'^ when i — >■ 0"*". Since 

when 77 — >■ first, it follows that f2{t) ~ for i — >■ 0+. Putting these 
together, assuming that f2{t) has the same critical exponent for both positive 
and negative t, yields 

7 = /3(5-l) (39) 
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which is a well-known relationship. To allow for the possibility of asymmetry 
about t = we write f2{t) = ±&± (±t)^, where the sign is chosen according 
to whether t is positive or negative, and the b± are both positive numbers.. 
The heat capacity at — )■ 



Cy ^ --/^t) 

c 

diverges as t^" when t — )■ 0^. Therefore we should write 



foit) = /o(t) - a+t 



2-a 



(40) 



(41) 



as t — ?■ 0+ where fo{t) is a smooth function. Hence the singular part of Cy is 
c+t^° with TcC+ = (2 — a)(l — a)a+. Another relationship among the critical 
exponents is 

a + 2/3 + 7 = 2 (42) 
Once again we allow for asymmetry about t = and write 



fo{t) 



/o(t)-a_(-t)2-" ift<0 
foit) - a+t^~" if t > 



(43) 



Using the best values arising from the 3D Ising model [23] and from data 
on real liquid-gas phase transitions one has /3 = 0.325 and 7 = 1.24. The 
above relationships then imply 6 = 4.815 and a = 0.11. Note that the mean 
field model considered previously respects both of these relationships among 
critical exponents too. 

Along the coexistence curve the chemical potential is 



(44) 



and the pressure is 



-2/2 (t) 



-fo{t) + fi{t) + \^Ut) 



-2/2 (t) 



(45) 



The formula for the isothermal compressibility along the coexistence curve 



IS 



(T-2 



-2(5-l)/2(t)(l±r;0' (46) 
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where the upper sign is for the hquid side and the lower sign is for the gas 
side {rjg = —rji). This indicates that /2 must of course be negative for t < 0. 
The formula for the heat capacity along the coexistence curve is 



Cy -- 

The singular part comes from the terms which are zero and second order in 
rj. This leads to cy — t- c_(— t)^" where 

2 

T, c_ = (2 - «)(1 - «)a_ + 7(7 - 1) (^^) ^- ■ (^S) 

According to Ref. [23] the thermal compressibility kt diverges as K^t~'^ 
when t — 7- 0"*" and as K^{—t)~'^ when t — )■ 0~, with ^ 5 a universal 

ratio. Also, the heat capacity at ?] — )■ diverges as c+t~" when t — )■ 0^ and 
as c_(— when t — >■ 0~, with Cj^/c_ ~ 0.5 another universal ratio. The 
former leads to the constraint 



(49) 



while the latter leads to 



2 

4-1 



2a+ = a_ + , , ( , 1 h- . (50) 

+ (2-«)(l-a) \afM) 

If we are not too far from the critical point we can use the following 
parameterizations. We can take /o- to be a constant. The function 

fu.\ m-b^{-tr ift<o 

m-^ f^{t) + h^v ift>o ^^^^ 

where f2{t) is a smooth function which vanishes at t = as a power bigger 
than 7. The function /i(t) is the chemical potential along the critical curve, 
which may be parameterized like this. Assume a quadratic relationship be- 
tween T and /i^,. 
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The curve hits the // axis at T = when ji — jiQ. The chemical potential at 
the critical point is 

/^c^/^oWl- (53) 



hence 



/i(t) = ne/.o^l-^(l + t)2. (54) 

It is apparent from the above expression for fi{t) that this whole parameter- 
ization is only good when T < Tq, otherwise fi{t) would become imaginary. 

Knowing this, it is straightforward to derive a simple formula for the 
latent heat per unit volume. Making use of phase equilibrium, AP = 0, 
T ^Tx, = i^x, Vg = -Vi, one finds 

Ae{t) = ^^^4t^ (55) 

which of course is valid only for t < 0. 

The independent constants may be taken as: a-,b-, fa, no, ^c, Tc, ric, plus 
the functions fo{t) and /2(t). This parameterization captures the critical 
exponents and amplitudes which are universal. To proceed further we require 
more information on the equation of state of QCD away from the critical 
point. 

For future reference it may be noted that one could add more terms to 
the free energy without changing the basic picture. For example, one could 

add f\{t)r]'^ and fQ{t)ri^. However, f4,{t) must vanish at t = so as not to 
affect the critical behavior. The coefficient fs{t) need not vanish at t = 
but it becomes irrelevant compared to the dominant term fa\ri\"' as rj 0, 
as do all powers of 77 greater than a. 



4 Fixing the Parameters 

In this section we narrow in on a phenomenological equation of state with in- 
put from various disparate sources. These include the results of lattice gauge 
theory calculations at zero baryon density, and models and extrapolations of 
the equation of state of cold dense nuclear matter. 

Concerning the function fo{t), what we know from the thermodynamic 
relations is that /o(0) = /Xc^^c — Pc = — TcSc and /o(0) = —TcSc- Hence for 
small t it starts out as fo{t) — — TcSc{l + t) + ---. 
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Suppose we wanted to extend this equation of state to T = 0. The 
minimum requirement is that the entropy density must vanish. This means 
that /^(-l) = /{(-I) = /^(-l) = /;(-!) = 0. This is satisfied by the 
above parameterizations of /i and fa, but /o and /2 must be modified. The 
simplest way to make /2(— 1) vanish which does not upset the critical point 
properties is to choose 

m = \h^it^ (56) 

following the idea that the function be smooth with the smallest integer 
power. Note that with this choice /2(i) is negative for all t < and positive 
for all t > 0. The simplest way to make /q(— 1) vanish, which does not upset 
the critical point properties, is to choose 

a_ = T,sJ{2 - a) . (57) 

If we then are so bold as to extrapolate our formula for the coexistence curve 
to T = we can input more physical information. 

For example, suppose we know that phase coexistence at T = occurs 
between a liquid phase with density ni{T = 0) and a gas phase with density 
ng{T = 0). Denoting normal nuclear density as no, we have no < ng{T = 
0) < ric < riiiT = 0) and ng{T = 0) + ni{T = 0) = 2nc. In that case we can 
solve for 6_ in terms of the density difference An = ni{T = 0) — ng[T = 0). 



a fa /An 



2-7 V2n, 



<5-l 



(58) 



Then fiQ = fi(T = 0) could be estimated by using the value calculated for 
nuclear matter at the density ng{T = 0). Now all the noncritical parameters 
would be determined apart from f^. 

Consider some common parameterizations of the cold nuclear matter 
equation of state [2l]. The energy density, pressure, and chemical potential 
are 

e = n{mN + Eo{n)) (59) 

P = n^^ (60) 
dn 

de , , 
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Case I: 



E^{n) 
H{n) 



K 
18 



n 



- 1 



no 



niN + -Eo(no) + 



K 



n 
no 



- 1 



n 



no 



(62) 



Case II: 

Eoin) 
fi{n) 



2K r 
~9" 



[n/no) 



1/2 



+ ^o(^o) 



niN + Eo{no) + 



2K 
~9~ 



[n/no) 



1/2 



2{n/noY^^ - ll (63) 



Here m^v = 939 MeV is the nucleon mass and Eo{no) = —16.3 MeV is the 
average binding energy per nucleon at nuclear matter density no = 0.153/fm^ 
[25] . The compressibility K is known to be 250 ± 30 MeV [251 126] ; we shall 
fix it at 250. See also [S]- 

Heavy ion collisions at the Bevalac and at the AGS showed no clear 
experimental evidence for the formation of quark-gluon plasma [28j. The 
baryon densities achieved were around two to four times nuclear matter den- 
sity. If one distributes one unit of baryon number within one electromagnetic 
radius of a proton, 0.8 fm, the baryon density would be about 0.47/fm^, 
which is slightly more than three times nuclear matter density. Therefore, it 
seems reasonable to estimate ng{T = 0) = 4no. At this density case I gives 
e(4no) = 641 MeV/fm^ and /i(4no) = 1381 MeV. Case II gives e{Ano) = 599 
MeV/fm^ and fi{4no) = 1089 MeV. Case I corresponds to a relatively stiff 
equation of state whose energy per baryon rises quadratically at high den- 
sity, whereas case II corresponds to a relatively soft equation of state whose 
energy per baryon rises linearly at high density. The energy densities are sim- 
ilar because they are dominated by the nucleon mass, not interactions. The 
chemical potentials differ by about 20% because interactions do contribute. 
The pressure is most sensitive to the interactions. Based on this information 
we estimate /io = 1230 ± 150 MeV. 

The Hagedorn temperature was already determined in the late 1960 's 
and early 1970's to be 160 MeV p9]. The critical temperature, no matter 
what order the transition is, ought to be slightly greater than this [30] . Data 
from heavy ion collisions at the SPS and RHIC show that no hadrons have 
ever been observed with a temperature greater than about 160 to 170 MeV 
(at very small chemical potential) [21]. Current lattice QCD calculations 
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agree that with the physical values of the light and strange quark masses, 
the transition is a rapid crossover at zero chemical potential. However, they 
disagree on the so-called critical temperature. One group [32] puts it at 150 
MeV while the other group [33] puts it at 190 MeV. Part of the discrepancy 
may be in exactly how this temperature is defined, but that is not entirely 
sufficient. Certainly, to accurately determine this temperature requires an 
accurate calculation of the low temperature hadronic equation of state. But 
this requires very fine lattice spacing, since the lattice must first discern the 
structure of individual hadrons, and it requires a very large lattice volume, 
since the hadrons become widely separated at low temperature. This is a 
difficult problem which may not be resolved for some time. However, it 
seems safe to estimate Tq = 170 ± 20 MeV. 

The value of the temperature at the critical point is of course not known. 
However, it should lie on or very near to the curve of T versus n under 
discussion. How then can we estimate the pressure Pc, energy density ec, 
entropy density Sc, and baryon density ric at the critical point? One obvious 
way is to use the formulas for a perfect massless gas of gluons and Nf flavors of 
quarks evaluated at Tc and fic- These formulas are {fi is the baryon chemical 
potential and quarks have one third of that value): 



e = 3P. (67) 



When the relationship 




/i^ = /x^ 1 - -I (68) 



with the aforementioned estimates of Tq and /io is used, it turns out that the 
pressure is almost independent of the numerical value of Tc. Since two phases 
in equilibrium with the same T and n have the same pressure, wouldn't it be 
nice if Pc was independent of Tc? This is one hint. A second hint is provided 
by the fact that all lattice QCD calculations show that the pressure, energy 
density, and entropy density are all lower than the ideal gas formula would 
suggest, at least at /i = 0. In fact, they indicate a negative contribution to 
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the pressure proportional to [31]-[3S], and a negative constant contribu- 
tion, like a bag constant. It has also been suggested, in the context of cold 
dense matter as might exist in neutron stars, that at T = there may be a 
contribution proportional to /i^ [36]. So let us hypothesize that, in the vicin- 
ity of the phase transition or crossover and above, the high energy density 
equation of state can be parameterized as 



P = A^T" + A2^^^T'^ + Aqu" - CT' - Djj,' - B 



(69) 



where A4, A2 and Aq are given by the perfect gas equation of state. Now 
suppose we substitute /i^ = — T'^/Tq) into this formula, and demand 
that it be independent of T. A simple exercise shows that Tq and fiQ must 
be related according to 



^ = 971^ 



i± A_ 

V 15 AhNf 



(70) 



Choosing Nf = 2.5, effectively to account for the smaller contribution from 
the heavier strange quarks at these temperatures and chemical potentials, 
and choosing the minus sign, leads to 



To 



6.67173. 



(71) 



in order that the coefficient of vanish. In other words, if To = 180 MeV 
then fiQ = 1209 MeV. This relationship is entirely in line with all the facts 
at hand. Demanding that the coefficient of the term vanish requires 



A2 - 2^Ao 



3.084T2 



(72) 



The lattice calculations of [33] found that 2C ~ 0.24 GeV^. This translates 
into C ~ 3.3Tq^ using their value of Tq ^ 190 MeV. There are no calculations 
of the fi"^ term in the pressure, but this analysis suggests that D is very small; 
we shall take it to be zero for simplicity of exposition. 

With even larger uncertainties Ref. [33] found that i? ~ Tq . For a rea- 
sonable interpolation of the lattice results near and just above the crossover 
region we take the coefficient to be 0.8. The parameterization is therefore 



„ 1697r2 5 a^a 

P = —rrr—T* + —^l^T^ + 



360 



36' 



6487r: 



3.084Tf?T2 



0.8To^ 



15 



1697r2 



+ —n^T - 6.168TJT 



s — 



90 



n 




(73) 



The critical pressure is computed from this to be Pc = 0.749Tq . The critical 
entropy density Sc, baryon density ric, and energy density ec of course depend 
on the choice of Tc and therefore Hc- 

All that remains is to specify An at T = and fa- Since fa is assumed 
to be constant, it is natural that it be proportional to Pc- For definiteness 
we shall take /o- = 5Pc ~ 512 MeV/fm'^ and An = nc/3. In what follows 
we shall use all of the above parameterizations and only vary Tc to see what 
effect it might have on heavy ion collisions. 



5 Numerical Results for Equation of State 



In this section we plot some of the thermodynamic functions that were de- 
rived in the previous two sections. It is important to note that, although 
the results do depend on the numerical values of the parameters, the critical 
behaviors obviously do not. In addition, since most of the parameters were 
chosen to match onto known properties of quark-gluon matter at /i — and 
to dense nuclear matter at T = the results should not be too far from what 
is at present impossibly difficult computations in QCD. 
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In Fig. 1 we show the pressure, entropy density, and energy density 
obtained from eq. (|75]) as functions of T at /i = 0. They are close to the 
curves computed in lattice QCD but not identical. The parameterization 
of eq. fl73l) is needed to extrapolate the lattice results to large chemical 
potentials. Anyway, all that is needed for the purposes of the chiral critical 
point are the values of Pc-, Sc and ec for a chosen value of T^, not the full T 
and fi dependence of the equation of state of quark-gluon plasma. 
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Figure 1: The pressure, entropy density, and energy density, normalized so 
that they all have the same asymptotic value, versus temperature at /i = 0. 
The parameterization is from eq. f l73l) which is motivated by lattice QCD 
calculations. 
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In Fig. 2 we show the curve of phase coexistence, T versus /x, according 
to eq. fl52p . One chooses somewhere along this curve. Then for T > Tc 
along this curve the transition is a rapid crossover, whereas for T < Tc along 
this curve the transition is first order. 

200 I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




1500 



/x(MeV) 

Figure 2: Temperature versus baryon chemical potential from the parame- 
terization of eq. (|52|) . The critical temperature lies somewhere along this 
curve. 
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In Fig. 3 we show the curve of phase coexistence, T versus n/no, for 
various choices of T^. These are indicated by the sohd curve. The dashed 
curve indicates the hmit of isothermal metastabihty, or isothermal spinodal. 
When T is scaled by Tc the phase coexistence curves fall on top of one 
another, as do the spinodals, indicative of a special scaling feature of this 
parameterization. 
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Figure 3: The solid curve denotes coexistence between high and low density 
phases. The dashed curve denotes the limits of metastabihty. When scaled 
by the critical temperature and density the curves lie on top of each other. 
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In Fig. 4 we show the isothermal compressibihty as a function of the 
temperature. For t < it depends on whether one approaches the phase 
coexistence curve from the low density side or the high density side. For 
i > it is computed at the critical density. The result is independent of the 
choice of Tc in this parameterization of the equation of state. 
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Figure 4: The isothermal compressibility. For t < they are evaluated along 
the coexistence curve while for t > it is evaluated at the critical density. 
The curves for = 60, 100, and 140 MeV lie on top of one another. 
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In Fig. 5 we show the heat capacity per unit volume as a function of 
temperature for various choices of T^. For t < it depends on whether one 
approaches the phase coexistence curve from the low density side or the high 
density side. For t > it is computed at the critical density. When scaled by 
the entropy density at the critical point, the result for i > is independent 
of the choice of Tc, whereas for i < it is almost but not quite independent. 
For reference, the entropy density at the critical point is 1.741, 3.416, and 
5.861 fm-3 for = 60, 100, and 140 MeV, respectively. 
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Figure 5: The heat capacity per unit volume. For t < they are evaluated 
along the coexistence curve while for t > it is evaluated at the critical 
density. When divided by the entropy density at the critical point the results 
are nearly independent of T^.. 
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In Fig. 6 we show the latent heat, or discontinuity in energy density, as 
a function of T for various values of T^. For 60 < < 140 MeV the latent 
heat is approximately 300 MeV/fm^ at T = and goes to zero at Tc. 
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Figure 6: The latent heat per unit volume versus temperature for three 
choices of Tc. 
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6 Fluctuations 



All thermodynamic functions are smooth and continuous for any finite vol- 
ume. The discontinuities associated with phase transitions only arise in the 
infinite volume limit. The question to be addressed here is whether the fea- 
tures characteristic of a chiral phase transition get smoothed out in heavy ion 
collisions to such an extent that they cannot be discerned from experimental 
observations. The analysis performed here is based on the Landau theory 
of fluctuations [37J, suitably modified to take into account the non- integer 
powers of rj appearing in the thermodynamic functions. 

In a uniform system of large but finite volume V, the thermodynamic 
potential Q depends only on the temperature and baryon chemical potential 
and is proportional to V. Due to finite size fluctuations at the given tem- 
perature and chemical potential, the actual value of the "order parameter" 
1] will not necessarily be the equilibrium one. To quantify this phenomenon 
we expand Q in powers of rj. 

T; 7]) = noifi, T) + r)r/ + fi2(/i, T)ri^ + T) |r/r (74) 

In equilibrium this must be an extremum with respect of variations in 77, 
namely 



^]l(^,T) + 2^]2(/i,T)r/ + afi,(/i,T)|r7r-^ign(r/) = 0. (75) 



The coefficient functions are determined by the fact that this condition is 
fulfilled by the equation of state. Noting the powers of rj which appear, it is 
clear that one should choose 



(76) 



where K is some factor yet to be fixed. Upon using the equilibrium relation 
between fi and rj, namely eq. (13T|) . one can determine the other coefficient 
functions. 



^^2 = Kf2 

= (77) 

Therefore 

n{fi,T-r]) = no{fi,T) + K\{-n,fi + f,)r] + f2v' + fa\vr] ■ (78) 
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In equilibrium Vt/V = —P. Comparing this with the expression fl3Up for the 
pressure, one can deduce that K = V and 



fio(/i,r) = r(/o-n,/i) . (79) 
Hence we have obtained the expansion around the equihbrium states 



n(/i, T; r]) = Qoifi, T) + V {-n,fi + f,)r] + f2rf + fM" 



(80) 



with flo{iJ,,T) given above. 

The probabihty V{ri) to find the system with a particular value of rj at 
given values of and T is 



~ e 



-n(M,T;r;)/T 



Along the coexistence curve n^/i = fi{t). Then 

n(/x, T- r,) - fio(/i, T) = V \hrf + fM" 



(81) 



(82) 



Since t < 0, f2{t) < and the thermodynamic potential has two equal min- 
ima at the densities of the liquid and gas phases, of course. This potential 
is shown in Fig. 7 for temperatures both below and above Tc. In this figure 
the volume was taken to be 400 fm^. Obviously the potential scales propor- 
tionately with this volume. The value of 400 fm^ is really quite optimistic 
for high energy nuclear collisions. Considering that the critical density is 
estimated to be about 5no ~ 0.75 baryons/fm^, this would mean that about 
300 baryons participate in the fiuctuation. That is a substantial fraction of 
the total of 394 in Au+AU, 416 in Pb+Pb, and 476 in U+U colhsions. Even 
then, the potential for the low and high density phases are only 5 MeV below 
the unstable mid-point when T/Tc = 0.6; it is even less as is approached. 

It is interesting to find the probability that the system has some value 
of rj other than rjg or t]i along the curve of phase coexistence. The relative 
probability is 



V{v)/V{r]i 



'Afl/T 



(83) 



where 



An = n{fi.,{T),T;r])-n{fi,{T),T-r]i) 

= v\fJv'~vf)+fA\vr-hn 



(84) 
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V 

Figure 7: The thermodynamic potential as a function of rj near the critical 
point; the volume is 400 fm^. The stable phases are located at the minima 
of the potential. Four different temperatures are shown, with the solid curve 
representing the critical temperature. 
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Figure 8: The probability to find the system at a particular density relative to 
the equilibrium densities at phase coexistence. Three different temperatures 
are shown. 

Still using V — 400 fm"', the relative probability is plotted as a function of 
f] in Fig. 8 for several values of T < Tc. For T > 0.6Tc there is more than 
a 90% probability to find the system with any value of 1] in the range from 
-0.2 to 0.2. A major reason that this probability distribution is so flat is due 
to the large value of the exponent a — 5.815 ^ 6. This is in contrast to the 
mean field models which have a = 4. For a smaller, probably more realistic 
volume from the perspective of nuclear collisions, the fiuctuations would be 
even greater. The magnitude of these fiuctuations suggests that it is difficult 
to probe the properties of the matter very close to the chiral critical point. 
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7 Conclusions 



In this paper we have constructed an equation of state in the vicinity of 
the chiral critical point. It incorporates the correct values of the critical 
exponents and amplitudes. Since only certain properties of the equation of 
state are universal, there is some freedom to vary the noncritical functional 
dependence on temperature and density and to change the parameters in 
those functions. The parameterization proposed here matches on to the 
equation of state at zero baryon density as calculated in lattice gauge theory, 
and at zero temperature using reasonable extrapolations of dense nuclear 
matter. Certainly, refinements and modifications are possible within the 
present framework. 

The Landau theory of fluctuations away from equilibrium states was 
employed to determine the potential magnitude of the fluctuations one might 
expect in heavy ion collisions. The magnitude of these fluctuations is quite 
large, partly due to finite volume effects but primarily because the critical 
exponent S is much larger than in standard mean field theories. This fiattens 
the Landau free energy as a function of density away from the equilibrium 
densities and hence decreases the cost to fluctuate away from them. 

In the future it would be highly desirable to have a parameterization of 
the equation of state that includes not only the behavior near the critical 
point but also extends to much higher temperatures and densities. Ulti- 
mately, to compare with experimental data, it will be necessary to incorpo- 
rate this knowledge into dynamical simulations of heavy ion collisions. 
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Appendix 



Here we review the parameterization of the equation of state near the chiral 
critical point as constructed in this manuscript for ease of apphcation. See 
the text for detailed explanations. 

The critical point lies somewhere along the curve 




= 1 



(85) 



Here Tq and Hq are constants. The pressure at the critical point is estimated 
from the expression 



TT 

90 



16 + 



2IN4 



+ 18^ + 3247r2 



^4 _Crp2_^ ^gg^ 



where Nf is the number of massless quark flavors. We use Nf = 2.5 to 
simulate the larger strange quark mass. The constants B and C are adjusted 
to represent the results of lattice QCD calculations in the vicinity of the 
crossover region at T > but /i — and to make the pressure a constant 
along the critical curve. Then 




(6.67173)^ 



(87) 



C 



32 



iVo /A 

18 V 15 45A?^/ 

B = o.8ro^ 



3.084T2 



(89) 

In particular Pc ^ 0.749rQ . The values of the entropy density, baryon den- 
sity, and energy density at the critical point are obtained from the above 
expression for the pressure via thermodynamic identities. When numerical 
values are required we use Tq = 180 MeV and thus /xq = 1209 MeV. 
The Helmholtz free energy is 



f-m+mv+f2{tw+fa{t)\vr 



(90) 



where r] = {n - nc)/nc and t ^ {T - T^)/T^. The value of a is 5.815. The 
coefficient functions are 




a_(-t)2-" ift<0 
a . if i > 



(91) 
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/i(t)=ne/.o^l-^(l + t)^ (92) 

fu^i h{t)-h.{-tr ift<o 

/o- = constant (94) 

The exponents are a — 0.11 and 7 = 1.24. The /o(t) and f2{t) are smooth 
functions of t. The critical amphtudes are related by 

K - (95) 

and 

Prom thermodynamic relations the smooth function 

Ut)^e,-T,s,{l + t) (97) 

to first order in t. The simplest parameterization of the other smooth function 
is 

h{t) = \h^^e (98) 

The parameters a_ and 6_ are 

a_ = T,s,/(2 - a) (99) 



(100) 

where An is the discontinuity in the baryon density at T = 0. For definiteness 
we use An = nc/3 and fa = 5Pc ~ 3.745Tq . 
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